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Periodic variations in the colours of the classical T Tauri star 
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Abstract. The classical T Tauri star RW Aur A is an irregular variable with a large amplitude in all photometric 
bands. In an extended series of photometric data we found small-amplitude periodic variations in the blue colours 
of the star, with a period of 2^64. The period was relatively stable over several years. The amplitude of the 
periodic signal is 0™21 in U-V, 0^07 in B-V, and about 0^02 in V-R and V-I. No periodicity was found in the 
V magnitude. The relevance of this photometric period to the recently discovered periodicity in spectral features 
of the star is discussed, and the hypothesis of a hot spot is critically considered. 
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1. Introduction 

RW Aur is a classical T Tauri star (CTTS) with an unusu- 
ally strong and variable emission line spectrum and large 
amplitude irregular variations in the brightness as well as 
in the veiling of the photospheric spectrum. The star is a 
resolved triple system (Ghez et al. 1993) with RW Aur A 
dominating the brightness in the optical spectrum. In a 
recent detailed investigation of RW Aur A (Petrov et al. 
|2001 , hereafter referred to as Paper I) periodic variations 
in many spectral features were found, with a period within 
2'?6-2'?9. The period was stable over several years of obser- 
vation. A model of non-axisymmetric magnetically chan- 
neled accretion seems to explain most of the observed vari- 
ations. The asymmetry of the accretion can be caused ei- 
ther by a close, low-mass, invisible companion (Gahm et 
al. 1999), or by a misalignment of the magnetic and rota- 



tional axes. 

The presence of an asymmetric accretion shock at the 
stellar surface should modulate the strengths of the Hei 
and He ii emission lines presumably formed in the post- 



shock region (Calvet & Gullbring 1998), while the hot pho- 
tospheric spot (10''^ K) below the shock should modulate 
the brightness of the star. In case of RW Aur A, clear pe- 
riodic modulations in Hei, Fell and other emissions were 
observed (Paper I) which gives a unique opportunity to 
study the expected periodic modulation of the brightness. 

In this paper we search for periodicities in the varia- 
tions of the brightness and the colours of RW Aur A, us- 
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ing all the available photometric data collected over three 
decades. 



2. Observational data 

The data were taken from the catalogue of U BVRI pho- 
tometry of TTS compiled by Herbst et al. ( |l994| ). For 
RW Aur A, the catalogue contains 575 observations in V 
and B-V, and fewer in other colours. The major part of the 
data are from the Majdanak Observatory (the ROTOR 
program by V.S.Shevchenko's group), and from Van Vleck 
Observatory (unp ublished) , other sources are Kardopolov 
& Filip'e v ( |1985[ ), Rydgren et al. (|l984| ) and Herbst et 
al. ( 1983| ). We also used the photometric data of 1996-99, 
presented in Paper I. Since the information about the ob- 
servational errors is not available in all cases, we adopt 
equal weights for all data. The light curve is shown in 
Fig. Ih, and the colour-magnitude diagrams are shown in 
Fig.0 

There is a good correlation between V and the colours 
B-V, V-R and R-I, with the slopes roughly correspond- 
ing to the mean law of interstellar extinction. It looks very 
much like the result of variable circumstellar extinction. 
Note, however, that for the spectral type of RW Aur A 
(K1-K4 V) the normal photospheric colour B-V should be 
0?"86-l"^05. With Ay = (Paper I) the observed pho- 
tospheric colour B-V is expected to be within CP95-1'.°14. 
Most of the observed B-V colours are much bluer. This ef- 
fect is also known for other CTTS: the stars typically have 
closest to normal photospheric colours when they are at 
minimum brightness, while becoming overly blue at max- 
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Fig. 1. Light curve of RW Aur A from 1965 to 1999 




The overall picture of the time variability of the star in 
the range of periods 0.8 — 1000 days is presented in Fig. |[ 
For the channels V, B and U, most of the variability is 
well concentrated around the zero frequency (the bunch 
of peaks around frequency 1 cyc/day is a result of aliasing 
from zero frequency). The relatively small peak at 2'?641 
in the U band is just the feature we are going to anal- 
yse in more detail. As it is well seen from the same figure, 
the strength of the peak for B-V is significantly amplified, 
and for the U-V colour (computed as (U-B)+(B-V)) the 
peak becomes the strongest in the full range of frequen- 
cies. Its statistical significance can be demonstrated using 
different methods. For instance, when we carried out 1000 
Monte-Carlo runs in the range of periods 2-3 days, the 
reshuffled data never showed peaks of this strength. For 
the full displayed range this gives a false alarm probability 
which is less than 0.0075. Due to the bad spacing of the 
data, any more detailed analysis can give only formally 
better levels of significance. This is why we tried to test 
the 2^641 feature by using a careful analysis of separate 
subsets of the full data. The peak at l'?609 in U-V (see 
Fig. H) is a result of aliasing (j-ggg « 1 — 9-545-), and we 
do not consider it further. 



Fig. 2. Colour-magnitude diagrams. The direction of red- 
dening is shown in the V versus B-V diagram according to 
the mean law of interstellar extinction. The UBVR bands 
are in the Johnson system, the I band is in the Cousins 
system 



imum light (Vrba et al. |1993|). This is believed to be due 
to an additional, variable non-photospheric hot continuum 
which also accounts for the veiling of the photospheric line 
spectrum. The continuum may belong either to hot pho- 
tospheric spots at the base of the accretion columns, or 
to plage-like regions of hot gas above the stellar surface. 
The contribution of the emission lines to the photomet- 
ric bands also accounts for a part of the excessively blue 
colours. 



3. Time series analysis 

3.1. Feature at P = 2} 6A1 

The compiled data set of RW Aur A is very inhomoge- 
neous. It contains large gaps and several dense subre- 
gions. To avoid any complication due to methodological 
artefacts, we selected for frequency analysis the simplest 
statistic available: the variance of the least squares fit 
residuals divided by the original variance of the data. As 
the model to be fitted we used a single harmonic curve 
/i + Acos{2TTtwk) + B sui{2Trt'Wk), where /i is the mean 
level of the light curve and Wk, {k — 1,2, ... , K) is one of 
the trial frequencies. The computed row of K variance ra- 
tios forms the least squares (LS) frequency spectrum (see 
Pelt |l992D . 



3.2. General and specific variability 

From the colour-magnitude diagrams (Fig. |^) we can con- 
clude that the total variability of RW Aur A consists of at 
least two components: general variability (GV), which is a 
large-scale variability in all the colours within about 0^"5, 
and specific variability (SV), which can be best seen in 
the colour-magnitude diagrams for U-B and B-V, as the 
scatter of the points around the linear trend. This scatter 
is very small in the red colours but steeply increasing to- 
wards the UV. Fig. || shows that the feature at P = 2^641 
belongs to SV (it is totally absent from the V band). In 
order to investigate the SV more closely, we should elim- 
inate the linear trend in the colour-magnitude diagrams 
and investigate the residual time variations in the colours, 
which can be defined as (B-V)ics=(B-V)-C-V-|-cowsi for 
B-V and in a similar way for other colours. Thus, we in- 
cluded an extra term into the LS fit of the harmonics: 
^ + Acos{2'Ktwk) + Bs\n[2TTtwk) + C ■ Y{t), and allowed 
coefficient C to vary together with the other model pa- 
rameters iJ.,A,B. As a result we obtained a new statistic 
for the period search which allows to optimally separate 
SV from GV. 

The results of such a procedure can be seen in Fig. ^ 
where short fragments of the LS spectra are depicted. One 
can see that a strong improvement is obtained for the B-V 
colour. The improvement is even more dramatic in the red 
colours, where the amplitude of the periodic signal is very 
small (SV<CGV). The phase diagrams with P = 2^641 for 
the residual variations (U-V)ros and (B-V)res are depicted 
in Fig. ^. The scatter is still quite significant so that we 
need to look at the data in more detail. 
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Fig. 3. Least squares spectra for magnitudes and colours. Due to the long time base, the actual spectra contain more 
than 300000 points each, so that the plot is somewhat schematic. For every pixel of the plot, the minimum and 
maximum value of the spectra covered by that pixel was computed, and the spectral content for the corresponding 
pixel was plotted as a vertical line from minimum to maximum. 
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Fig. 4. Details of the LS spectra around period P — 2^641. 
Specific values of the test statistic (denoted here by D) are 
also displayed. It is well seen that the peaks are amplified 
in the spectra of the residual variances. 
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Fig. 5. The phase diagram for (U-V)ros and (B-V)i.es with 
the period 2^641. The best fit harmonic is also given as a 
continuous line. 



3.3. Stability of the period 

Here, we will focus on the analysis of the residual variation 
(U-V)res- The largest number of observations are concen- 
trated in the following four seasons: 1981-82, 1987-88, 
1988-89 and 1989-90. In Fig. | we plotted the results of 
the LS fits of a single harmonic with the period P = 2'?641 
for each of the four seasons separately. The scatter for two 
seasons (1981-82 and 1988-89) is reasonably small, but for 
season 1987-88 it is much larger, and for the last season 



with the small number of points it is very high. However, 
the amplitudes (maximum-minimum) are still significant 
(0^33, a'^20, 0^40, a"15, respectively), especially if com- 
pared to the amplitude 0™21 of the full data fit (see Fig. ||). 
There is also a reasonably good phase stability: the phase 
shifts around the full fit phase are 20%, 5.7%, 3.3%, 8.3%, 
respectively. 

As can be seen from Fig. ^ the feature around the 
period 2^641 consists of two main depressions in the LS 
frequency spectra. The smaller peak right from the main 
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, Local fits for (U-V)ics with the period 2.641 days, 
are 59 points for the 1981-82 season, 54 points for 
!8, 54 points for 1988-89, and 27 points for 1989-90. 



peak (P = 2'?639) is a result of a cycle count ambigu- 
ity between seasons 1981-82 and 1987-88. If we com- 
puted the LS spectrum for the three subsets from 1981- 
1989 together, ignoring all other data points, the peaks 
around 2^641 and 2^639 occured to be comparably deep. 
Consequently, the last digit in the value of the accepted 
period depends on wildly scattered points outside the 
above mentioned seasons and must be taken only as the 
best approximation. 

In conclusion, we can say that the most prominent 
feature in the frequency spectrum of the U-V colour is 
certainly a real regularity which is more or less persistent 
at least during eight years (see Fig. ^. 

4. Discussion: Is there a hot spot? 

The periodic variations of the blue colours are most prob- 
ably related to the periodicities in the spectral lines dis- 
cussed in Paper I. The spectroscopic data set used in 
Paper I covered only four years (1996-1999), which al- 
lowed to find a period within the range 2'?6 — 2^9. The 
period of 2^77 was selected as the best one for demon- 
stration of the phase diagrams in Paper I, although with 
a period of 2^64 the diagrams look similar. In this paper, 
using the much larger time span of the photometric data 
we confirm the period and prove that it was relatively sta- 
ble over at least several years. The stability of the period 
probably excludes the possibility that the periodicity is 
caused by a local asymmetry in the magnetospheric struc- 
ture, which is supposed to be a short-lived phenomenon. 
We still cannot distinguish between the binary hypothesis 
and the hypothesis of the inclined magnetic rotator (mis- 
alignment of the magnetic and rotational axes) , discussed 
in Paper I. 

The amplitudes of the periodic colour variations ris- 
ing steeply towards the UV indicate the presense of a hot 
source of radiation (hot spot). When the whole data set 
is used, the full amplitudes of the sinusoidal oscillations 
are: 0"^21 in U-V, 0™07 in B-V, and about (r02 in V-R 
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Fig. 7. Phase diagrams with the period 2'?641. Upper: flux 
(in 10~^^ ergcm"^ s~^) in the narrow Hel emission. The 
errors are due to uncertainties in the decomposition of 
the He i emission profile into the narrow and broad com- 
ponents. Centre: equivalent width of the Nai Di accre- 
tion component. Lower: the residual B-V colour (flUed cir- 
cles: the photometry made simultaneously with the spec- 
troscopy; open circles: the whole set of the photometric 
data). 



and V-I. Using black body approximations for the fluxes 
from the photosphere and the hot spot, one can estimate 
the temperature of the hot spot and the fraction of the 
visible stellar disc covered by the spot (Vrba et al. 1993). 
With Tphot=4800K, and the amplitudes given above, we 
get Tspot=15000K and the covering fraction 0.1%. Such 
a small spot gives only a 0™05 amplitude in the V magni- 
tude. 

In Paper I we presented strong arguments in favour 
of the magnetospheric accretion model. The accelerated 
gas stream(s) flowing toward the star is clearly seen in 
the red-shifted components of many spectral lines. The 
strengths of the red-shifted components vary periodically, 
being increased considerably at the phase when the ac- 
cretion stream is in the line of sight. We found that the 
narrow He I emission varies in correlation with the red- 
shifted accretion components, and concluded that the nar- 
row emission originates from the shock region(s) at the 
bottom of the accretion stream, near the stellar surface. 
The broad emissions are believed to be formed in the 
global magnetosphere of the star. The two relevant phase 
diagrams are reproduced here with the new period 2^641 
and compared to the colour phase diagram (Fig. 

We might now identify the small hot spot with the 
accretion shock at the stellar surface. Then we would ex- 
pect the bluer colours when the narrow He i emission is 
stronger. However, from Fig. one can notice that quite 
the opposite is observed: when the colours are bluer, the 
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Fig. 8. Fluxes in the narrow and broad components of the 
He I emission vs. the residual B-V colours. 



flux in the narrow He I emission is lower. This inverse cor- 
relation is better seen in the Fig. || with the logarithmic 
scale for the flux. The " normal" correlation exists between 
the colour and the broad line flux: the larger the broad line 
flux, the bluer the colour (see also Paper I). 

Hence, the colour varies in " anti-phase" to the strength 
of the accretion components and to the flux in the narrow 
He I emission, in the sense that the star is redder when the 
shock region is facing the observer. There is no doubt that 
magnetospheric accretion is going on in RW Aur A, and 
that the accretion column(s) exists, but the expected hot 
spot at the base of the accretion column does not reveal 
itself in our series of data. The small hot spot, discussed 
above, may belong to the non-uniformely distributed hot 
gas, which also radiates in the broad emission lines. 

The apparent absence of the hot spot related to the 
accretion column is really puzzling. Either the accretion 
shock is absent indeed, or the complicated geometry of 
the flows and possible extinction of light within the gas 
streams makes the observed variations so complicated. 
Simultaneous spectroscopic, photometric and polarimet- 
ric monitoring of the star could clarify the situation. 
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